Insulating Behavior of a Trapped Ideal Fermi Gas 
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We investigate theoretically and experimentally the center-of-mass motion of an ideal Fermi gas in 
a combined periodic and harmonic potential. We find a crossover from a conducting to an insulating 
regime as the Fermi energy moves from the first Bloch band into the bandgap of the lattice. The 
conducting regime is characterized by an oscillation of the cloud about the potential minimum, while 
in the insulating case the center of mass remains on one side of the potential. 

PACS numbers: PACS: 03.75.Ss Fermi gas, 03.75.Lm BEC in optical lattices, 63.20.Pw localized states, 
05. 45. -a Non linear dynamics 



Recent experiments have demonstrated the possibility 
of reaching the superfluid regime with trapped fermionic 
atoms 0,0- In this context, the combination of trapped 
cold fcrmions and optical lattices is very promising. 
On the one hand it represents the natural analogy with 
superconductivity of electrons in crystal lattices, while on 
the other it might be a additional tool to manipulate the 
interaction properties of the system and to investi- 
gate the superfluid phase. Atomic Fermi gases in lattices 
are interesting also in view of observing novel quantum 
phases that have been predicted for interacting Fermi pj 
and Fermi-Bose systems j^J- 

In this work we investigate a Fermi gas of non interact- 
ing atoms trapped in a harmonic field and subjected to a 
one-dimensional optical lattice. We study both theoreti- 
cally and experimentally the center-of-mass dynamics of 
the system after a sudden displacement of its equilibrium 
position along the lattice. We find that the statistical dis- 
tribution makes the dynamical properties of such system 
highly non trivial. The most dramatic new feature is that 
the system can exhibit a conducting or an insulating be- 
havior, depending on whether the Fermi energy lies in 
the first Bloch band of the lattice or in the gap region. 
In the first case the cloud oscillates symmetrically in the 
harmonic potential, while in the second case it remains 
trapped on one side of the potential. By tuning the width 
of the lattice band we observe the full crossover between 
the two different regimes. 

Semiclassical Model. Let us first introduce our theo- 
retical model to describe the system. Neglecting the in- 
teratomic collisions, the many-body Hamiltonian can be 
simply written as a sum of single particle Hamiltonians: 
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where s is a tunable dimensionsless parameter, Er — 
3=r(^f) is the recoil energy and A is the wavelength of 
the laser creating the lattice. Since the harmonic oscilla- 
tor length is typically much larger than the lattice spac- 
ing d = A/2 we can study the system in the semiclassical 
approximation. 

We first concentrate on the z-component of the Hamil- 
tonian In the limit of a vanishing harmonic 
field (uj z — > 0), the eigenfunctions of the stationary 
Schrodinger equation are Bloch waves. The cigenener- 
gies s n (p z ) have the typical Bloch band structure as a 
function of the quasimomentum p z , defined in the first 
Brillouin zone (—irh/d < p z < +nh/d), and of band 
index n (in the following we consider only the lowest en- 
ergy Bloch band n = 1). In a semiclassical approach, we 
can incorporate the effects of the periodic potential mod- 
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ifying the kinetic part of the Hamiltonian: — > e(p z ), 
while the harmonic confinement generates a driving field. 
Therefore, Eq. Q can be approximated as: 

H = e( Pz ) + \mLo 2 z z 2 + + X - mU j\{x 2 + y 2 ). (2) 

The effects of a sudden small displacement Zq of the 
center of the harmonic trap can be included adding the 
perturbative term H pert (z, t) = miu 2 Q(t)zoz, where 9(t) 
is the unit step function. A first insight in the behav- 
ior of the system is provided by the study of the linear 
dynamics in a 1-D configuration. 

1-D Analysis. The physics becomes particularly clear 
in the semiclassical phase space, as shown in Fig. ^ The 
lines are isoenergetic single particle orbits; they belong 
to two different classes, separated by the dashed orbit, 
which has an energy 2S, corresponding to the width of 
the first Bloch band. The orbits are close when their 
energy is in the band, i.e. the particles oscillate around 
the trap minimum. The orbits are instead open when 
their energy is in the gap, because the quasimomentum 
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FIG. 1: Phase trajectories for a trapped 1-D Fermi gas in a 
lattice at T — 0, just before and after the displacement of 
the trap (figs. A and B, respectively), and their dynamical 
evolution (figs. C and D). The ordinate and abscissa are in 
units of p z = p z d/2h and z = \J muJlz 1 /AS. 



is reversed at the band edges. Particles in these orbits os- 
cillate on the sides of the harmonic potential, performing 
the analogous of Bloch oscillations in a linear potential. 
In general the oscillation frequency of both kinds of orbits 
shows a strong dispersion with energy. 

As shown in Fig.^^j the Fermi gas at T — uniformly 
fills the phase space region with energy below Ep. A sud- 
den displacement of the center of the harmonic potential 
corresponds to a shift of the center of the phase space (see 
Fig.njJ). The blue region contains particles that are still 
in equilibrium in the new configuration of the trapping 
field. The red and yellow regions, containing particles 
moving on close and open orbits respectively, are instead 
out of equilibrium and give rise to a collective dipole mo- 
tion. The phase space region opens and melts during the 
dynamics as a consequence of the energy dependence of 
the single particle oscillation frequency (see figs. C, D), 
yet leaving constant the phase space volume (because of 
the Liouville theorem) and therefore preserving the Pauli 
principle. To one hand, the red orbits dephase on a longer 
time scale with respect to the yellow ones. Therefore, 
the relaxation and the frequency of the oscillation mode 
are dominated by the particles moving around the cen- 
ter of the phase space, in the red region. To the other 
hand, the yellow orbits remain open, trapping the center 
of mass of the system on one side of the harmonic poten- 
tial. Notice that in this 1-D configuration, the damping 
of the oscillation disappears in the linear limit (small ini- 
tial displacement) at T=0. 

3-D Analysis. The semiclassical phase-space distribu- 
tion f(r,p,t) is governed by the Liouville equation 
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scribes the quantum dynamics up to orders of 0(h 2 ) 0. 
For small displacements, we can linearize the distribution 
function: f(r,p,t) = f (f,p) + g{r,p,t), where f (r,p) = 

^exp (f3(Ha(r,p) — /i(T))) + lj is the fcrmionic solu- 
tion of the stationary Liouville equation for temperature 
T relative to Hq, and g(r,p, t) is a small time-dependent 
correction induced by the perturbation H pert {z,t). The 
dipole oscillations evolve according to : 
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where Q(x(u;)) is the imaginary part of the dipole re- 
sponse function x( w )- I n order to calculate it we solve 
Eq. J2Jl in the linear regime by the general method illus- 
trated in , by using the energy of motion in z-direction 
E z = e(p z ) + ^moj z z 2 as an independent variable instead 
of p z . The imaginary part of the linear response function 
then becomes 
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Such an equation correctly de- 
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In the above equations u>o(E z ) — loq(E z ,z — Z2{E Z )) is 
the single particle frequency, zi(E z ) and zi{E z ~) are the 
inversion points of the motion (i.e. the maximum and the 
minimum positions reached by the particle of energy E z ). 
Notice that, due to the 3-D nature of the system, all val- 
ues E z < Ep contribute to the response function (J5J even 
at zero temperature, where the integral inside the square 
brackets becomes (E z — Ep)Q(Ep — E z ). As a conse- 
quence a damping due to dephasing is in general present 
also in the linear limit, which is different from the 1-D 
case. As shown by the Eqs. (14171 . in order to calculate 
the dipole motion we need to know explicitly the Bloch 
energy e{p z ) as a function of the quasimomentum. In 
tight binding approximation (large lattice height, corre- 
sponding to small tunnelling between the potential bar- 
riers) the first band is e(p z ) — 25sm 2 (p z \/4K) where 
8 = 5(s) [2( is a decreasing function of the lattice height. 
With this choice for the dispersion the frequency ojq(E z ) 
can be expressed in terms of elliptic integrals [lfj and we 



3 



0,5 - 



-0.5 - 



; ; 

1 


i > 




; • 


Q 


> 9 ' ' 


> ; t 


> J i 


J 1 Y« 




v 4/ ' • 0-^. 





_l Lz j ' j ' , ^ ^ 

50 100 150 200 250 

Time (ms) 

FIG. 2: Dipole oscillations of the Fermi gas of 40 K atoms 
at T = 0.3 Tf in presence (filled circles and full line) and in 
absence (empty circles and dotted line) of a lattice with height 
s = 3. The lines are the theoretical predictions, the circles 
are the experimental results. The horizontal dot-dashed line 
represents the trap minimum. 

can calculate analytically the quantities of interest. In 
regimes out of tight binding we use a parametrization of 
the energy band [llj. 

Experiment. We realize experimentally the system by 
using a Fermi gas of 40 K atoms. The sample is prepared 
in the F = 9/2, nip = 9/2 state and sympathetically 
cooled |l2j to quantum degeneracy in a magnetic trap 
with frequencies ui z — 2tt x 24 s _1 and ojj_ = 2tt x 275 
s . Subsequently a one-dimensional optical lattice is 
adiabatically superimposed along the weak axis of the 
trap (z direction) The lattice is created by a laser 
beam in standing-wave configuration, with wavelength 
A = 863 ran. The lattice height can be adjusted in the 
range U = 0.1 - 8E Rl where E R /k B = 317 nK. The typ- 
ical Fermi gas is composed by 25 000 atoms at a temper- 
ature that can be varied between 0.2 Tp and Tp, where 
the Fermi temperature is Tp w 300 nK. The system is 
brought out of equilibrium by displacing the magnetic 
trap minimum along the lattice. The typical displace- 
ment is Zq=15 /itm, much smaller than the 1/e 2 radius 
of the cloud (110 /xm) 13]. After a variable evolution 
time in the trap the atoms are released from the com- 
bined potential. We detect the position of the center of 
mass of the cloud by absorption imaging after a ballistic 
expansion of 8 ms. 

Results. In Fig. @ we show both the theoretical 
prediction (solid line) and the experimental observation 
(solid circles) of the dipole oscillation of a Fermi gas at 
T = 0.3 Tp in presence of a lattice with s=3. For com- 
parison, we show also the dipole motion in the absence 
of the periodic potential (dashed line and open circles) 
which consists in an undamped harmonic oscillation at 
the trap frequency u) z . We observe a dramatic change 
of the dipole oscillation in presence of the lattice, which 
is well described by the 3-D model. The appearance of 
an offset in the oscillations is due to the significant frac- 
tion of particles moving along open orbits; for the given 
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FIG. 3: Relative oscillation Relative offset of the oscillations 
of a Fermi gas in the lattice (normalized to the initial displace- 
ment) as a function of the lattice height. The circles are the 
experimental data and the continuous line is the theoretical 
prediction for 3xl0 4 atoms at T=100 nK. The dashed line is 
the prediction for a Fermi gas of 2.5xl0 3 atoms at T=0. 



parameters the Fermi energy is indeed larger than the 
bandwidth 25 w OAEp. This fraction behaves macro- 
scopically as an insulator, because its center of mass does 
not move under the harmonic force but stays trapped on 
one side. The fraction of the gas occupying closed or- 
bits can instead oscillate in the harmonic potential, and 
has therefore a conducting nature. A damping appears 
as expected because of the dephasing between different 
orbits. Also, the oscillation frequency is reduced because 
of the larger effective mass of the atoms in the lattice. 
To understand how such phenomenology depends on the 
width of the first energy band, we have performed a se- 
ries of measurements by keeping the atom number and 
temperature of the Fermi gas constant, and varying the 
lattice height. In this way it was possible to move the 
Fermi energy within the energy gap between the first 
and second band of the lattice. In Figs. 13141 we plot the 
measured dependencies of the offset, damping rate and 
oscillation frequency, that are in good agreement with 
the theoretical calculations. Fig. [3] shows the crossover 
from a conducting behavior in low lattices to an almost 
completely insulating behavior in higher lattices. Here 
we have plotted the relative oscillation offset defined as 
z osc / ' zq, where z osc is the center of oscillation of the sys- 
tem in the lattice, as a function of the lattice height. 
Note how the relative offset, which represents the insu- 
lating fraction of the Fermi gas, stays small as long as 
2<5 < Ep, and then raises quite rapidly towards unity. 
Since in the present experiment Ep « E R , an insulat- 
ing fraction appears already with low lattices; the the- 
ory however shows that the threshold for the insulation 
moves to higher lattices in case of smaller Fermi energies 
(dashed line in Fig. EJ) . The disagreement between ex- 
periment and theory at low lattice heights, s < 3, can 
arise from the population of higher bands due to the fi- 
nite temperature and/or Landau-Zener tunnelling, which 
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are not included in the theoretical model. Increasing s, 
the energy gap between bands increases, and the single 
band calculations become more realistic. 

In Fig. 0| we show the observed features of the con- 
ducting fraction of the gas. The damping rate of the os- 
cillation shown in Fig. also increases with the lattice 
height, because of an increased dispersion of the oscilla- 
tion frequencies of atoms in closed orbits As shown 
in Fig. 0Jd, the oscillation frequency of this component 
of the gas is close to that expected for a particle at the 
bottom of the band u>o(E z = 0). Actually, the theo- 
retical analysis shows that this value is strictly reached 
only asymptotically during the oscillation, because ini- 
tially closed orbits with different frequencies contribute 
to the dipole motion. 

As a consequence of the Pauli principle, which keeps 
the energy distribution broad, the Fermi gas exhibits an 
insulating behavior even at T=0. The observed phenom- 
ena only weakly depend on the gas temperature, at least 
in the region 0.2-1 Tp that we have explored so far in 
both experiment and theory [hi ], and in general we ob- 
serve an increase of both offset and damping for increas- 
ing temperatures, as expected because of the broader en- 
ergy distribution. An analogous insulating phenomenon 
can be observed in uncondensed Bose gases at tempera- 
tures larger than the bandwidth. Indeed also in this case 
a significative part of the single particle orbits will have 
energy in the gap regi on g iving rise to insulation. How- 
ever we have observed [l5| that in the case of bosons, or 
even in the case of fermions admixed with bosons, inter- 
atomic collisions quench the open orbits and eventually 
bring the system into the equilibrium position. 

Conclusion. In this paper we have investigated the 
dynamics of a non-interacting Fermi gas in presence of 
a combined periodic and harmonic potential, comparing 
the experimental observation with a semiclassical theory. 
We have shown that the system crosses from a conduct- 
ing to an insulating behavior by tuning the Fermi energy 
relatively to the first energy band of the lattice. The 
present work provides a basis for further investigation 
of interacting Fermi gases in lattices, in both normal and 
superfluid phases, including the superfluid Josephson-like 
oscillations already observed in a Bose-Einstein conden- 
sate confined in the same trapping potential 16]. 
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FIG. 4: A) Comparison between theory (line) and experi- 
ments (circles) for the damping rate of the dipole oscillations 
of the Fermi g function of the lattice height. B) Oscil- 

lation frequency of the Fermi g function of the lattice 

height. The line is the expectation for a particle oscillating 
at the band bottom. 
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